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Recent progress in the determination of both masses and radii of neutron stars are starting 
to place stringent constraints on the dense matter equation of state. In particular, new theoretical 
developments together with improved statistical tools seem to favor stellar radii that are significantly 
smaller than those predicted by models using purely nucleonic equations of state. Given that the 
underlying equation of state must also account for the observation of 2Mq neutron stars, theoretical 
approaches to the study of the dense matter equation of state are facing serious challenges. In 
response to this challenge, we compute in a model-independent way the underlying equation of 
state associated with an assumed mass-radius template similar to the “common radius” assumption 
used in recent studies. Once such a mass-radius template is adopted, the equation of state follows 
directly from the implementation of Lindblom’s algorithm; assumptions on the nature or composition 
of the dense stellar core are not required. By analyzing mass-radius profiles with a maximum mass 
consistent with observation and common radii in the 8 to 11km range, a lower limit on the stellar 
radius of a IAMq neutron star of Rns ^ 10.7 km is obtained in order to prevent the equation of 
state from violating causality. 

PACS numbers: 26.60.-c, 26.60.Kp, 21.60.Jz 


How does subatomic matter organize itself and what 
phenomena emerge is one of the overarching questions 
guiding the field of nuclear physics [T] . In the case of 
atomic nuclei, the quest to answer this question requires 
understanding the nature of the nuclear force and the 
limits of nuclear existence. In the case of extended nucle¬ 
onic matter, this involves elucidating the nature of neu¬ 
tron stars and dense nuclear matter. In this letter we 
focus on the latter. 

Owing to the long-range nature of the Coulomb force, 
extended nucleonic matter must be electrically neutral. 
As a result, dense nuclear matter must be by necessity 
neutron-rich. This is because the electronic contribution 
to the energy increases rapidly with density, so electron 
capture becomes energetically advantageous. Given that 
such extreme conditions of density and isospin asymme¬ 
try can not be realized in terrestrial experiments, neutron 
stars have become unique laboratories for the exploration 
of dense matter. This situation has created a strong syn¬ 
ergy between nuclear physics and astrophysics, that has 
been cemented even further through an intimate inter¬ 
play between theory, experiment, and observation [2]. In¬ 
deed, powerful telescopes operating at a variety of wave¬ 
lengths drive new theoretical and experimental efforts 
which in turn suggest new observations. 

A recent example of such a unique synergy is how 
accurate measurements of massive neutron stars El S] 
have informed nuclear models that fall under the gen¬ 
eral rubric of density functional theory. Density func¬ 
tional theory (DFT) offers a comprehensive—and likely 
unique—framework to describe strongly interacting nu- 
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clear many-body systems ranging from finite nuclei to 
neutron stars. Rooted on the seminal work by Kohn 
and collaborators [5] , DFT shifts the focus from the com¬ 
plicated many-body wave function to the much simpler 
one-body density. The implementation of DFT to nuclear 
physics requires that the parameters of the model—which 
encode some of the complicated many-body dynamics— 
be determined by fitting directly to experimental data. 
In this regard, the accurate measurement of neutron star 
masses has been vital to the accurate calibration of some 
modern energy density functionals EHS]. 

However, whereas the determination of neutron star 
masses is accurate and beyond question, attempts to re¬ 
liably extract stellar radii EHII] have been hindered by 
large systematic uncertainties that resulted in an enor¬ 
mous disparity in stellar radii—ranging from as low as 
8 km E] all the way to 14 km m- It appears, however, 
that since those first analyzes were performed, the sit¬ 
uation has significantly improved in the last few years 
through a better understanding of systematic uncertain¬ 
ties, important theoretical developments, and the im¬ 
plementation of robust statistical methods [MIS]. Al¬ 
though a consensus has yet to be reached, these re¬ 
cent studies seem to favor stellar radii in the 9-11 km 
range. Particularly intriguing among these are the re¬ 
sults by Guillot and collaborators that suggest a “com¬ 
mon radius” of i?NS = 9.llj;5km for all five quies¬ 
cent low mass x-ray binaries used in their analysis m-. 
this common-radius value has been slightly revised to 
Rns = (9.4± 1.2) kmdSj. What makes this result espe¬ 
cially provocative is that satisfying the small radius and 
large mass constraints simultaneously is enormously chal¬ 
lenging. Indeed, to our knowledge no optimized energy 
density functional can simultaneously reproduce both of 
these constraints. And from the very large number of 
models available in the literature we are aware of 
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only a few that account for both [TBl - ED] . From these, 
the model due to Wiringa, Fiks, and Fabrocini [18] re¬ 
lies on a microscopic approach based on the Argonne U 14 
(AV14) nucleon-nucleon potential supplemented with the 
Urbana VII (UVII) three-nucleon potential. This model 
predicts a maximum stellar mass of 2.13 Mq and a “com¬ 
mon radius” of i?NS —10-4km. The other theoretical ap¬ 
proach due to Hebeler and collaborators mnoi is also 
microscopic in nature, but instead uses nuclear interac¬ 
tions derived from chiral effective field theory. In partic¬ 
ular, their softest equation of state is consistent with the 
^2Mq limit and predicts a radius as low as iii .4 = 9.7km 
for a I. 4 M 0 neutron star [2(1] . 

In this letter we aim to elucidate the nature of neutron 
star matter by relying on a powerful result first proven 
by Lindblom more than two decades ago [21] . It is a well 
known fact that all stellar profiles may be determined 
from the Tolman-Oppenheimer-Volkoff (TOV) equations 
once an equation of state P = P{£) relating the pressure 
P to the energy density £ is supplied. That is, given an 
equation of state (EOS), the TOV equations 
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can be solved once a value for the central pressure P{0) = 
Pc and enclosed mass M( 0 ) = 0 are specified at the origin. 
In particular, the stellar radius R is determined from the 
condition P{R) = 0 and the corresponding stellar mass 
as M = M{R). In this manner, the EOS generates the 
mass-radius (MR) relationship for neutron stars. What 
Lindblom was able to prove and implement is that the 
inverse also holds true: knowledge of the MR relation 
uniquely determines the neutron star matter equation of 
state [ 2 T] . 

In the spirit of Lindblom’s approach, we now proceed 
to outline the methodology required to obtain the unique 
equation of state associated with an assumed MR tem¬ 
plate that closely resembles the “common-radius” hy¬ 
pothesis adopted in Ref. [H] . This choice of template is 
the only model assumption made in this work; no addi¬ 
tional assumptions on the nature or composition of dense 
matter are required. For alternative attempts at obtain¬ 
ing generic model-independent constraints see Ref. [22] 
and references contained therein. 

The particular MR template adopted in this work fol¬ 
lows from our latest optimization: a relativistic energy 
density functional labeled “FSUGarnet” [H]. Such a tem¬ 
plate, among many others depicted in Fig[T| is displayed 
with the purple solid line containing a few additional cir¬ 
cles. FSUGarnet was optimized with inputs from the 
ground-state properties of finite nuclei, their monopole 
response, and a maximum neutron star mass consistent 
with observation. Yet, even after such an optimization, 
extrapolations to the high-density domain remain highly 


uncertain. However, it is worth underscoring that rela¬ 
tivistic density functionals provide a Lorentz covariant— 
and hence causal—framework that becomes critical as 
one extrapolates to the high densities encountered in 
the stellar interior. Of particular relevance to this work 
are predictions for the maximum neutron star mass and 
radius of a 1.4Mq neutron star; FSUGarnet predicts: 
Mmax = (2.07 ± 0.02) Mq and R 1.4 = (13.0 ± 0.1) km. 
Although such a stellar radius appears to significantly 
exceed some of the preferred limits, the predicted EOS is 
relatively soft in the immediate vicinity of nuclear matter 
saturation density. This is relevant as the pressure in the 
neighborhood of twice nuclear matter saturation density 
sets the overall scale for stellar radii [ 13 - 



fig. 1: (Color online) Mass-Radius profiles as predicted by 
the microscopic model of Ref. [IS] and by the four relativistic 
density functionals considered in the text. As a check, the 
circles on top of the FSUGarnet profile were obtained using 
an EOS extracted from Lindblom’s algorithm. Also shown 
are the four profiles R8-R11 generated from the FSUGar¬ 
net template and having radii of 8-11 km, respectively. The 
horizontal band shows the observational constraints from the 
mass measurements reported in Refs. Oil], whereas the ver¬ 
tical band displays the constraint on stellar radii published in 
m- Also shown are regions excluded on purely theoretical 
grounds: BH is the black hole limit, L&P is the causality con¬ 
straint obtained in Ref. m, and the grey region defines the 
causality constraint obtained in this work. 

Having defined the underlying template, additional 
MR curves may be easily generated by simply shifting the 
full FSUGarnet profile to the desired radius. Indeed, the 
curves labeled R8-R11 in Fig. [^ were generated in pre¬ 
cisely this manner and are characterized by stellar radii 
of 8-11 km (rather than 13 km) for a 1 . 4 M 0 neutron star. 
In the case of these four curves, the underlying EOS is 
unknown. Besides the FSUGarnet template, MR predic¬ 
tions generated with other relativistic density function¬ 
als are also included for comparison. These range from 
NL3 [231 [23] and FSUGold [2S] that were optimized with¬ 
out incorporating neutron star information, to the more 
recent FSUGold2 parametrization that includes the 2Mq 
constraint in the calibration [7] . Note that all these func¬ 
tionals provide an accurate descriptions of a variety of 
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ground-state properties of finite nuclei. Finally, displayed 
with the red line is the MR relation predicted by the mi¬ 
croscopic formalism of Wiringa, Fiks, and Fabrocini [TB]. 
With the exception of FSUGold, all MR relations are con¬ 
sistent with the 2Mq limit depicted in the figure with the 
narrow magenta band. However, none of the relativistic 
density functionals satisfy the common-radius constraint 
of Ref. [TS] depicted in Fig.j^by the vertical (cyan) band. 

Our next step is to use Lindblom’s approach to ob¬ 
tain the EOS corresponding to the R8-R11 profiles dis¬ 
played in Fig.[^ Briefly, the “inversion” procedure is 
implemented as follows. First, one assumes complete 
knowledge of the EOS up to a pressure Pi and energy 
density £i capable of supporting neutron stars up to a 
mass of about 0.4 M©. Specifically, using such an EOS 
and a central pressure of Pc = Pi generates, through the 
TOV equations, a neutron star of mass Mi^OA Mq and 
radius Ri. Second, one steps along the mass-radius tra¬ 
jectory by selecting a neutron star of mass > Mi 

and radius i?i+i whose central pressure Pc and energy 
density £c are to be determined. To compute Pc and 
£c, the TOV equations are integrated inwards from the 
surface to the core, with boundary conditions given by 
P{Ri+i) — Q and M(Ri_|_i) = M^+i. Hydrostatic equilib¬ 
rium guarantees that as one moves towards the interior 
of the star, the pressure will continue to increase until the 
maximum known pressure Pi will be reached at a certain 
radius Ti located near the center of the star. Given that 
Ti is close to the origin, one may use suitable series ex¬ 
pansions to determine the central pressure Pc and energy 
density £c in terms of known quantities, i.e., mass, pres¬ 
sure, and energy density, at the small radius [^. In 
this manner, after one iteration the EOS is extended from 
{Pi,£i) to (Pi+i,£i+i) = {Pc,£c)- In the next iteration, 
one proceeds in exactly the same fashion, namely, one 
moves another step along the mass-radius trajectory and 
repeats the algorithm with the newly augmented equa¬ 
tion of state. After the whole mass-radius trajectory is 
sampled, the entire high-density component of the equa¬ 
tion of state is mapped out. Testing the reliability and 
accuracy of Lindblom’s inversion algorithm is fortunately 
very simple. To start, one uses a known EOS—such as 
the one predicted by FSUGarnet—to compute its asso¬ 
ciated MR relation. Then, one applies Lindblom’s algo¬ 
rithm to extract the “new” EOS directly from the MR 
profile. Finally, one verifies that both the new EOS and 
the resulting MR profile are consistent with the originals. 
The result of such a test is illustrated in EigJ^for the case 
of ESUGarnet; the circles represent the results obtained 
with an EOS extracted from Lindblom’s algorithm. 

Having verified the accuracy of Lindblom’s algorithm, 
we can now obtain the equations of state that generate 
the R8-R11 MR profiles assumed in Eig.[2 The resul¬ 
tant equations of state along with those predicted by the 
non-relativistic and relativistic models are displayed in 
Fig.d Note that among the relativistic density function¬ 
als, FSUGarnet displays the softest EOS at low energy 
densities (f < 200 MeV/fm^). Eor reference, FSUGarnet 



FIG. 2: (Color online) Equations of state associated to the 
mass-radius profiles displayed in Fig.[^ In the case of the 
R8-R11 profiles, the equations of state were obtained from 
Lindblom’s algorithm m- Also shown is the P = 8 line which 
demarcates the boundary beyond which an equation of state 
becomes ultrabaric. Equations of state that cross the line 
become superluminal at considerably lower energy densities. 


predicts a pressure for pure neutron matter at saturation 
density of Ernm = (2.60±0.08) MeV/fm^. Given that the 
pressure in the vicinity of twice nuclear matter saturation 
density sets the overall scale for stellar radii [T7], FSUG¬ 
arnet generates the mass-radius profile with the smallest 
radii; in contrast, NL3 with the stiffest EOS generates the 
largest stellar radii. However, since the maximum neu¬ 
tron star mass is sensitive to the EOS at higher densities, 
ESUGold (with Mmax = 1.72 Mq) becomes the softest be¬ 
yond £~200MeV/fm3. 

Also displayed in Fig.is the P — £ line which marks 
the boundary beyond which an equation of state becomes 
ultrabaric Eauz]. Being rooted on a Lorentz covariant 
framework, the predictions from all relativistic density 
functionals lie safely below the ultrabaric line. However, 
this is not the case for the R8-R11 profiles. Could this be 
an indication that such profiles violate causality? That 
is, could the speed of sound exceed the speed of light in 
the high-density cores of these extremely compact stars? 
To test if a given EOS respects causality, we compute 
the speed of sound in the medium c^/c^ = dP/d£ and 
determine whether a point exists at which the EOS be¬ 
comes superluminal. Such a point {Pc,£c) determines 
the central pressure and energy density of the heaviest 
neutron star that can be supported by the given EOS; 
beyond such a mass, the EOS required to support the 
star becomes acausal. 

For clarity, we display in Fig.|^ the various equations 
of state on a reduced linear scale; only NL3 is included as 
a representative member of the relativistic density func¬ 
tionals. The figure clearly indicates that with the excep¬ 
tion of NL3, all five equations of state become superlu¬ 
minal much before they cross into the ultrabaric region. 
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FIG. 3: (Color online) As in Fig.|^ but now on a reduced lin¬ 
ear scale, we display some of the equations of state considered 
in this work; for simplicity, only NL3 is included as a repre¬ 
sentative member of the relativistic density functionals. With 
the exception of NL3, which is causal at all densities, all other 
equations of state become superluminal at the pressure and 
energy density indicated by the blue squares. In turn, these 
values determine the central pressure and energy density of 
the heaviest star that may be supported by a causal EOS. 


The pressure and energy density at which the EOS be¬ 
comes superluminal are depicted by the blue squares. In 
turn, these values determine the central pressure and en¬ 
ergy density of the most massive neutron star that may 
be supported by a causal EOS; beyond this value the 
mass-radius relation becomes unphysical. We found that 
the value of this critical mass Merit decreases rapidly as 
the common radius decreases. Specifically, although in 
the case of the Rll template we find Merit = 2.06 Mq, 
and thus consistent with the present 2Mq limit 1313 , 
we obtain Merit/MQ = 1.83,1.57,1.26, for RIO, R9, and 
R8, respectively. Note that the microscopic model of 
Wiringa, Fiks, and Fabrocini violates causality beyond 
Merit « 1.96 Mq [H] . These values of Merit define the 
lower boundary of the grey region labeled as “Causality 
this work” in Fig.[2 We thus conclude, within the scope 
of the adopted MR template, that in order to satisfy the 
current 2M0 limit, the stellar radius of a \AMq neu¬ 
tron star must exceed 10.7 km. This represents the main 
finding of our work. 

To place our newly developed constraint on the proper 
context, we display in Fig.[l] two other commonly used 
limits. The weakest of the two (labeled “BH”) represents 
the black hole limit that precludes a neutron star of mass 
M from having a radius R equal to its Schwarzschild ra¬ 
dius. This constraint has no impact on the mass-radius 
profiles considered in this work. A more stringent con¬ 
straint (labeled “L&P”) follows from an analysis by Lat- 
timer and Prakash m that uses a realistic EOS up to a 
certain value of the energy density that is then matched 
to a causally limited EOS. By doing so, they obtain a 
limit on the stellar compactness that may be written 


as R > 2.83GM/c^. Although stronger than the BH 
constraint, this causality limit does not affect the re¬ 
cent analysis by Guillot and Rutledge m- Indeed, the 
L&P limit for a 9 km neutron star is 2.16 M©. On the 
other hand, our analysis provides an even more strin¬ 
gent constraint of 1.57 Mq, and pushes the minimum ra¬ 
dius to 10.7 km—barely consistent with the upper limit 
quoted in Ref. |15j , but fitting comfortably within the 
10.8('1 q' 4 range given in the very recent work by Ozel 
et al. [16] . 



FIG. 4: (Golor online) Sensitivity to small variations to the 
standard MR template. Here the RIO template (dashed line) 
is modified to produce two additional MR profiles, both with 
a radius of 10 km for a IAMq neutron star. The right portion 
of the plot displays the critical mass beyond which the under¬ 
lying EOS becomes superluminal. Finally, the labels indicate 
the minimum radius required to support a 2Mq neutron star 
with a causal EOS of the given shape. 

Given that the MR profile is the only assumption made 
in this work, it is pertinent to ask how sensitive are 
our conclusions to the underlying shape. As shown in 
Fig.ll a constant-radius profile is characterized by a large 
(indeed infinite) derivative in the mass region of about 
0.6 ;< M/Mq <1.6. If this shape is changed by keep¬ 
ing the slope large but now negative (i.e., an increase 
in mass is accompanied by a reduction in the radius) 
then the compactness will increase and the underlying 
EOS will become superluminal even earlier. However, if 
instead the slope becomes positive (i.e., an increase in 
mass is followed by an increase in radius) then it be¬ 
comes harder to assess whether such an MR profile will 
generate an acausal EOS; such a shape is displayed by 
the green solid line in Fig.|^ Still, our results suggest 
that if variations to the standard template are not overly 
dramatic, our main conclusions remain valid. Indeed, the 
minimum stellar radius got shifted by only O.lkm; from 
10.7 to 10.6 km. 

In summary, motivated by the latest progress in mea¬ 
suring neutron star radii, which combined with measure¬ 
ments of massive neutron stars provide critical insights 
into the dense matter equation of state, we have exam- 
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ined whether neutron stars may be as compact as re¬ 
cently suggested. To do so, we relied on a powerful al¬ 
gorithm developed by Lindblom to obtain the equation 
of state from knowledge of the mass-radius relationship. 
That is, given a MR profile, the underlying EOS may 
be obtained without any assumption on the nature or 
composition of the dense stellar core. Thus, all model 
dependence lies in the assumed MR profile, which in the 
present study was chosen to have a nearly constant ra¬ 
dius shape. Using such a template, MR profiles were 
constructed with a maximum mass of 2.1 Mq and with 
common radii spanning the 8 to 11 km interval. Using 
Lindblom’s algorithm, the equation of state associated 
with each of these MR profiles was obtained. Further, by 
imposing causality, namely, by enforcing that the speed of 
sound be less than the speed of light, we obtained a strin¬ 


gent constraint on the maximum compactness of neutron 
stars. Indeed, by demanding that a causal EOS be able 
to support a 2Mq neutron star, we obtained — within 
the scope of the adopted template — a lower limit on the 
stellar radius of a 1.4M0 neutron star of 1?““ = 10.7 km. 
Note that our result imposes a lower limit on stellar radii. 
In contrast, recent observational studies are placing up¬ 
per limits on neutron star radii. We trust that such a 
theoretical-observational synergy will continue to prove 
beneficial in the coming years in order to determine the 
dense matter equation of state. 

This material is based upon work supported by the 
U.S. Department of Energy Office of Science, Office 
of Nuclear Physics under Award Number DE-FD05- 
92ER40750. 
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